Transient Regime of Kerr Frequency Comb Formation 
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Temporal growth of an optical Kerr frequency comb generated in a microresonator is studied both experi- 
mentally and numerically. We find that the comb emerges from vacuum fluctuations of the electromagnetic field 
on timescales significantly exceeding the ringdown time of the resonator modes. The frequency harmonics of 
the comb spread starting from the optically pumped mode if the microresonator is characterized with anoma- 
lous group velocity dispersion. The harmonics have different growth rates resulting from sequential four-wave 
mixing process that explains intrinsic modelocking of the comb. 

PACS numbers: 42.62.Eh, 42.65.Hw, 42.65.Ky, 42.65.Sf 



Kerr combs excited in nonlinear optical microresonators 
hold promise as chip scale generators of octave spanning op- 
tical frequency combs with unique characteristics lU. They 
also reveal a rich, and as yet not well understood, variety of 
nonlinear dynamical phenomena. Different regimes of comb 
generation have been observed and several theoretical mod- 
els developed. While the mainfocus of the research in these 
studies is related to understanding the spectral properties of 
the Kerr comb, less work has been devoted to its time do- 
main behavior. Several experimental and theoretical |j^ 
1^ studies of mode locking regimes of this nonlinear process 
have been published very recently, but the problem of tem- 
poral growth of comb harmonics was addressed only theoret- 
ically 10, In this Letter we report on experimental study 
of transient dynamics of Kerr frequency combs and also pro- 
pose a theoretical explanation for the observed results. We 
also suggest a numerical model, backed by the experiment, 
that predicts a different dynamics for Kerr comb growth com- 
pared with earlier theoretical predictions. Our research clearly 
explains why harmonics of the comb are modelocked. 

Kerr combs are generated from electromagnetic vacuum 
fluctuations due to modulation instability of a continuous 
wave (cw) light confined in an externally pumped nonlinear 
dispersive resonator When the power of the cw optical pump 
that is nearly resonant with one of the modes of the resonator 
exceeds a certain threshold, the cw field inside the resonator 
becomes unstable, and multiple frequency harmonics are gen- 
erated in the modes. The harmonics are equally spaced due 
to energy and photon number conservation laws, imposed by 
four-wave mixin g pr ocess (FWM), which is responsible for 
comb generation Oiol [Till . 

The growth of the comb is not instantaneous. It was found 
10, S that formation of a fully developed comb can take up 
to a hundred ring-down periods of the resonator mode. The 
goal of the present contribution is to measure this time inter- 
val directly. We performed a numerical simulation and found 
that the DC power of light exiting the resonator depends on 
the degree of comb formation. The effective intrinsic qual- 
ity factor of the pumped optical mode suddenly drops long 
after the steady state amplitude of the circulating light in the 
mode is reached. The reduction of the intrinsic quality fac- 
tor changes the attenuation of the pump mode since the bal- 
ance between intrinsic loss and coupling loss in the resonator 
changes. Hence, the time delay between the moment the light 



enters the pump mode and the moment when the frequency 
comb is generated can be directly measured by detecting the 
light escaping the resonator It is also possible to track the 
time dependence of the power of the RF signal generated by 
the comb to measure the delay. We performed such a mea- 
surement with a Kerr frequency comb produced in a high-Q 
CaF2 whispering gallery mode (WGM) resonator, and experi- 
mentally confirmed the theoretical prediction of |0,[^ as well 
as the result of our numerical simulations. In what follows we 
describe our simulations and the experiment in detail. 

To simulate the transient regime of the Kerr comb we use 
the theoretical model developed in llsl [T2ll . We introduce an 
interaction Hamiltonian 

V = ~^{e^fe\ (1) 

where 



is the coupling parameter obtained under assumption of com- 
plete space overlap of the resonator modes, wo is the value of 
the optical frequency of the externally pumped mode, c is the 
speed of light in vacuum, n2 is the cubic nonlinearity of the 
material, V is the effective geometrical volume occupied by 
the optical modes in the resonator, and no is the linear index 
of refraction of the resonator host material. The operator e is 
given by the sum of annihilation operators of the electromag- 
netic field for 41 interacting resonator modes that we took into 
consideration: 

41 

The external cw pump is applied to the central mode of the 
group, so the Kerr frequency comb is expected to have twenty 
red- and twenty blue-detuned harmonics with respect to the 
frequency of the pumped mode. We take into account only 
the second order frequency dispersion that is recalculated for 
the frequency of the modes. 

Equations describing the evolution of the field in the res- 
onator modes are generated using Hamiltonian (flji and input- 
output formalism developed for ring resonators lUSll 

o-o = -(70 + i^j)aj + J[t>, hj] + i^oe"*'^*<52i,j, (4) 
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where 821 j is the Kronecker's delta; 70 = 70c + 7oi is the half 
width at the half maximum for the optical modes, assumed to 
be the same for the all modes involved; and 70c and 7oi stand 
for coupling and intrinsic loss. The external optical pumping 
is given by 



hujn 



(5) 



where P is the value of the cw pump light. We neglect the 
quantum effects and do not take into account corresponding 
Langevin noise terms. 

The set (|4| should be supplied with an equation describing 
the light leaving the resonator. Assuming the pump power 
does not depend on time, the relative amplitude of the output 
field is given by ifTsIl 



^out 



270c 7oe 
7oc + 7oi ^0 ' 



(6) 



where tq = 27ri?no /c is the light round trip time for the res- 
onator, and R is the radius of the resonator. It is assumed that 

e{t - To) w e{t). 
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FIG. 1: (a) Optical frequency comb generated in the WGM res- 
onator for selected numerical parameters of the system, (b) Time 
dependence of the normalized amplitudes of the first five harmon- 
ics of the optical frequency comb generated in the WGM resonator 
characterized with (.g/70)^''^ = 5 x 10~^. 

We solved Eqs. (|4]i and (|6]l numerically taking into account 
the interaction of 41 optical modes. Since the number of 
modes in our numerical simulation is limited, we selected a 
resonator with large GVD to be sure that the spectral bound- 
ary conditions do not impact the nonlinear process. We as- 
sumed 2^21 - W22 - W20 = -70, (-Po/7o)(g/7o)^/^ = 1-5, 
and Lo = UJ21 — 1.770, and found that the resonator generates 
the optical frequency comb with spectrum shown in Fig. ([T^). 
The comb has a frequency harmonic in each second optical 
mode (for the selected values of the parameters of the sys- 
tem), which sometimes is observed experimentally (see Fig. 5 



in fl^). The different regimes of Kerr comb generation will 
be studied elsewhere, while in this work we focus on the tran- 
sient processes. 

We performed a numerical simulation of the growth of the 
first five harmonics of the frequency comb Fig. ([TJ)) revealing 
several important features of comb generation. The sidebands 
grow exponentially with different growth rates. The first har- 
monic ((1) in Fig. [U has the slowest growth rate of 70, the 
second harmonic ((2) in Fig. [T]l has twice faster growth rate, 
270, the third harmonic - 870. The first harmonic starts to 
grow much earlier than others, the second harmonic starts be- 
fore the third, etc. Such a temporal behavior shows that the 
particular realization of the Kerr comb is initiated by hyper- 
parametric oscillation 1 12] that involves only two optical side- 
bands closest in frequency to the pump. The next order of 
optical harmonics is generated in the stimulated comb due to 
FWM of the already generated sidebands and the pump light 
Ih']. This stimulated process does not have a threshold. In 
other words, the threshold of comb generation coincides with 
the threshold of the hyper-parametric oscillation. The lowest 
order harmonics also determine the phase as well as the fre- 
quency of the rest of the harmonics. Such a comb is always 
phase locked and optical pulses are formed in the resonator 
(Fig.EJ. 
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FIG. 2: Relative power and phase for the optical pulses leaving the 
resonator. The shape of these pulses is similar to the shape of the 
pulses reported in | 2]. 

Another important observation is related to the temporal be- 
havior of the pump light. Since the growth of the comb har- 
monics is not instantaneous, the pump light is not influenced 
by the comb growth initially. The pump power is impacted 
only after the frequency harmonics approach their saturation 
values. This behavior is clearly seen at the phase diagram 
(Fig.[3t), where the pump mode has two attractors. The first 
one ((I) in Fig. [3^) corresponds to the steady state solution for 
for the pump light in the nonlinear resonator with no harmon- 
ics generated, and the second attractor ((II) in Fig. [3^) corre- 
sponds to the steady state solution with the saturated comb. 

The duration of the transition process depends on the non- 
linearity parameter 5/70 which defines the maximal number 
of photons generated in the comb harmonics. The larger g is, 
the smaller is the number On the other hand, the initial pho- 
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ton number is equal to unity. The growth rate of the comb is 
on the order of 70. Therefore, the transient process is longer 
when g /70 is smaller. We have calculated the delay for two 
cases: (5/70)^^^ = 5 x 10~^ (estimated for a small resonator, 
e-g- l Uil) and (5/70)^''^ = 10^^ (estimated for a much larger 
resonator). The result of the calculation is shown in Fig. ^p). 

The behavior of the DC power of light leaving the resonator 
has a certain peculiarity. The value of the power exiting the 
resonator initially decreases and then increases (Fig.lSji). The 
phenomenon can be explained from the stand point of criti- 
cal coupling ifTsll . There is no light at the output of a linear 
resonator if 70^ = 70c and the steady state is reached; all the 
pump light is absorbed in the resonator That is why the ampli- 
tude of the exiting light drops after the pump is on. A number 
of the pump photons is redistributed between harmonics of the 
comb, as the comb is generated. Those harmonics leave the 
resonator reducing the absorption of the pumping light. The 
interference phenomenon resulting in the critical coupling is 
also deteriorated since the pump light confined in the corre- 
sponding mode changes its amplitude and phase due to the 
nonlinear process. 
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FIG. 3: (a) Transient behavior of the normalized amplitude of 
the electromagnetic field of the externally pumped optical mode, 
(021 (t)). Since the mode is initially empty, its amplitude increases 
and reaches the attractor (I) approximately during resonator's ring 
down time, after the cw pump light is on. Then, after a certain 
period of time equal to the time interval needed for the comb to 
grow, the amplitude drops to a certain level and changes its phase 
approaching another attractor, (II). (b) The transient behavior of 
the light exiting the resonator is calculated for slightly overcoupled 
modes, 70c ~ 1.57oi, and for different values of nonlinearity g. 
The oscillator reaches its steady state faster for larger nonlinearity 
(id/lo)^^^ ~ 5 X 10^^, solid blue line), and slower for smaller 
nonlinearity (((?/7o)^'^^ = 10~*, dashed red line). 

The temporal dependencies shown in Fig. ([T})) and Fig. 
can be verified experimentally if one measures the power of 
the comb harmonics exiting the resonator Instead of the direct 
measurement of the optical harmonics, the power of the radio 



frequency (RF) signal generated on a fast photodiode by the 
comb can be measured. We performed such experiments. 

We used a calcium fluoride (CaF2) whispering gallery 
mode (WGM) microresonator with loaded Q-factor 2.2 x 10^ 
(full width at the half maximum of the mode is 90 kHz, and 
corresponding ring down time 1.75 /is). The intrinsic Q-factor 
of the resonator was 5.5 x 10^ (70c — 1.57oi), which means 
that the attenuation of light in the modes was primarily given 
by the interaction with the evanescent field coupler (a glass 
prism), and not by the scattering and loss of the resonator host 
material. The resonator was pumped using a 1545 nm dis- 
tributed feedback (DFB) semiconductor laser, self-injection 
locked to a selected resonator mode 111 61 . The laser was emit- 
ting approximately 6 mW of power, 30% of which entered the 
resonator. 
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FIG. 4: Schematic of the experimental setup. The WGM resonator 
is pumped with a DFB laser, self injection locked to the selected 
mode. The output light is analyzed with an optical spectrum ana- 
lyzer (OSA) showing the spectrum of the generated Kerr comb. Part 
of the light is sent to a fast photodiode (PD). The photocurrent, mod- 
ulated with a frequency equal to the comb repetition rate, is directed 
to an RF power detector (RFD), and a fast oscilloscope. This sig- 
nal is proportional to the convolution of the harmonics of the optical 
frequency comb. Some of the light is also detected with a slow pho- 
todiode (PD) and the photocurrent from this photodiode is forwarded 
to another channel of the same oscilloscope. This signal shows the 
integral DC power leaving the resonator. Comparing the signals with 
the fast oscilloscope we are able to measure the time delay between 
the generation of the Kerr comb and the moment the pump light en- 
ters the corresponding mode. 



We tuned the laser frequency by changing the laser current, 
and observed the temporal behavior of the signals at the OSA 
and the oscilloscope (Fig. |4|. We noticed a significant de- 
lay between the start of the optical pumping of the mode, and 
the moment of comb generation (Fig. |5]). Some experimen- 
tal data are very similar in shape to the theoretical predictions 
(Fig.|5^,b), while other data are not (Fig.|4j;,d). However the 
delay is always there. The different shapes of the signals are 
measured since we selected different modes in the resonator. 
It was shown that the GVD of the modes change depending 
on the resonator morphology ifTHl and depending on the ex- 
ternally pumped mode. In addition to the GVD, the coupling 
efficiency is also different for different modes. 
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FIG. 5: Experimentally observed transient behavior of the Kerr fre- 
quency comb measured via monitoring the power of the DC signal 
on a slow photodiode generated by the light exiting the resonator 
(curves (a) and (c)) as well as the power of the RF signal generated 
by the comb on a fast photodiode (curves (b) and (d)). Two modes 
are considered. The coupling to the mode resulting in the comb that 
produces in curves (a) and (b) is much lower compared with the cou- 
pling to the mode responsible for curves (c) and (d). The modes are 
also characterized with different GVD values. 

It may appear that conditions of the numerical simulation 
and the experiment described above are not exactly the same. 
In the experiment we used a large CaF2 resonator (6.72 mm 
in diameter) with free spectral range of 10 GHz. The GVD 
of modes of an ideal spheroidal resonator of this size is much 
smaller compared with the value of GVD we utilized in the 
numerical simulations. Moreover, GVD is normal for the fun- 
damental mode sequence of the CaF2 resonator (assuming its 
ideal spheroidal shape). Nevertheless, the optical frequency 
comb observed in the experiment has a spectral shape similar 



to the one obtained with the simulation, and the transient pro- 
cesses measured in the experiment have good correspondence 
with the theory. This contradiction is resolved if we take into 
account the notion that the morphology of a monolithic res- 
onator allows changing the sign and value of GVD [11]. In 
the experiment we did not use the ideal spheroidal resonator 
or the fundamental mode sequence. We rather sent the light 
to a higher-order mode of the resonator. The observed Kerr 
frequency comb was generated for light tuned nearly at the 
top of the WGM resonance (a; ~ luq). This is possible only 
if GVD is anomalous and large for a locally selected mode 
family, in accordance with the theory of hyper-parametric os- 
cillation 1 12]. That is why the numerical simulation for a res- 
onator characterized with large anomalous dispersion is appli- 
cable for the description of the experiment performed with a 
resonator seemingly having normal dispersion. 

To conclude, we have studied the transient regimes of Kerr 
frequency comb formation in a nonlinear monolithic optical 
resonator. We found that the well developed comb genera- 
tion is delayed by tens of ring down intervals of the resonator. 
Kerr combs created in larger and/or less nonlinear resonators 
have longer transient period compared with those generated in 
smaller and/or more nonlinear resonators. We noted that Kerr 
combs generated in resonators possessing comparably large 
anomalous GVD are always mode-locked. Finally, we ex- 
perimentally found that there exist mode families with large 
anomalous GVD even in WGM resonators made out of a ma- 
terial with normal GVD. We experimentally validated the re- 
sults of our numerical simulations with a suitable mode fam- 
ily. 
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